Coulomb-induced Rashba spin-orbit coupling in semiconductor quantum wells 
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In the absence of an external field, the Rashba spin-orbit interaction (SOI) in a two-dimensional 
electron gas in a semiconductor quantum well arises entirely from the screened electrostatic potential 
of ionized donors. We adjust the wave functions of a quantum well so that electrons occupying the 
first (lowest) subband conserve their spin projection along the growth axis (s z ), while the electrons 
occupying the second subband precess due to Rashba SOI. Such a specially designed quantum well 
may be used as a spin relaxation trigger: electrons conserve s z when the applied voltage (or current) 
is lower than a certain threshold V"; higher voltage switches on the Dyakonov-Perel spin relaxation. 
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Introduction. The ability to control the amplitude of 
the Rashba[l] spin-orbit interaction (SOI) electrically 
[3|] makes this type of SOI one of the most promising in- 
struments for manipulating spins of electrons in future 
spintronic devices[4, 5]. The most commonly studied 
case is the Rashba SOI in a two-dimensional electron gas 
(2DEG) inside a symmetrical quantum well. In such a 
system the Rashba SOI is zero in the absence of an exter- 
nal electric field, applied perpendicularly to the plane of 
the quantum well. Further, in a quantum well with more 
than one subband, such external-field induced spin-orbit 
coupling does not show strong dependence on the sub- 
band index. 

It is natural to ask whether Rashba SOI can also appear 
in the absence of an external electric field. A well-known 
argument, [6j based on the Ehrenfest theorem, leads to 
the conclusion that this cannot happen, as long as the 
confinement potential for the valence band is propor- 
tional to the confinement potential for the conduction 
band (which is usually the case). However, this argu- 
ment fails to take into account the screened electrostatic 
potential of the donors: this breaks the proportionality 
between the potentials acting on electrons and holes and 
leads to a finite SOI in the absence of an external field. 
We describe this effect as Coulomb-induced Rashba SOI. 

An attractive feature of the Coulomb-induced Rashba 
SOI is that it can be made strongly subband-dependent 
by a proper engineering of the shape of the quantum 
well. In this letter we use techniques of inverse scatter- 
ing theory[7] to find a shape of the quantum well for 
which the first (lowest) subband is free of SOI, while 
the second has a rather large SOI. We suggest that such 
a shape may be realized by digital alloying techniques. 
In such a specially engineered quantum well, electrons 
in the lowest subband would conserve their spin pro- 
jection s 2 along the growth axis z ([110]), but electrons 
in the second subband would suffer spin relaxation by 
Dyakonov-Perel mechanism. More importantly, we sug- 
gest that a parallel electric field (as opposed to the conven- 
tional perpendicular electric field) might be used to control 




FIG. 1: The proposed device. The 2DEG is confined inside 
the quantum well grown along the 2- ([110]) direction; the 
current flows along the x-direction. The widths of the spacers 
separating the central zone (the well) from the doping layers 
on the two sides of the well are different. 

the SOI amplitude by "pumping" electrons into the sec- 
ond subband. The electrostatic potential generated by 
this non-equilibrium population would "turn-on" the 
Rashba coupling, leading to Dyakonov-Perel spin relax- 
ation also in the lowest subband. 

The model. We consider a two-dimensional elec- 
tron gas confined within a quantum well (QW) in an 
Al T Gai_ v As heterostructure (see Fig. [I). 

We assume that the electrons in the QW are provided 
by donors outside the well, and that the concentration of 
donors is high enough so that the 2DEG is degenerate. 
The electronic states in the well are plane waves with 
momentum p in the plane of the well (x, y). They are 
further characterized by a subband index n, which de- 
termines the transverse envelope wave function (p n (z). 
Our QW contains only two subbands for electrons, n = 1 
or 2, and only the lowest subband (n = 1) is populated 
in equilibrium. 

If the quantum well is grown along the [110] direction, 
the Dresselhaus[8] spin-orbit interaction is absent, and 
the main source of spin relaxation is the Rashba[9] spin- 
orbit interaction, 

H R = (X n (0 X py - Oyp X ), (1) 
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where the SOI amplitude a n is given by 



a„ = (n\d z \lv\n) 
3n 
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U v is the completef^ potential acting on the holes|0], 
Eg = 1 .52 eV is the fundamental band gap, A = 0.34 eV, is 
the spin-orbit splitting of the valence band, P = 10.5 eVA 
is the matrix element of momentum between s and p 
atomic states, and \n) are subband envelope functions. 
Apart from the bare valence band confinement E v and 
a possible external potential U ex t (such as the one that 
arises from a perpendicular electric field), U v includes 
the self-consistent "Hartree" potential Uh induced by 
the donors and the inhomogeneous charge distribution 
inside the well: 



U V (Z) = !iext + E v (z) + lT H (z). 
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For x < 0.45, the confinement potential in the valence 
band of Al x Gai_. v As is proportionally to the confine- 
ment potential, E c (z), in the conduction band: 

E v (z) = 2£ c (z)/3. (4) 

Setting Uext - and making use of Eq. <(4j) we can write 

U v = ^U c (z) + ^!i H (z), U c (z) = E c (z) - LTh(z), (5) 

where U c (z) is the full potential acting on the electrons. 
The expectation value of the force <9 2 HJ C in the state \n) 
vanishes by the Ehrenfest theorem, so that only the con- 
tribution from the Hartree potential survives. We con- 
clude that, in the absence of external fields, the Rashba 
SOI is produced entirely by the Coulomb contribution Uh 
to the hole potential, 
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where U\\ is obtained from the self -consistent solution of 
both Schrodinger and Poisson equations for the electron 
envelope functions (p n (z). It is evident from ((6) that the 
SOI vanishes in all subbands if the confinement potential 
E c (z) and the donor distribution po(z) are symmetric. 
The reason is that in this case d z Uu(z) is an odd function 
of z, while the envelope function (p n (z) has a definite 
parity. 

Spin-orbit coupling engineering. When the confine- 
ment potential E c and/or the donor distribution po is 
asymmetric, both SOI amplitudes a\ and a 2 are, in gen- 
eral, non-zero. In order to eliminate the Dyakonov-Perel 
relaxation in equilibrium, we have to engineer a confine- 
ment E c for whichjl2|l ct\ = but a 2 + 0. This can not be 
achieved by breaking the symmetry of either E c or po', 
one has to break both symmetries, so that 



E c (-z) + E c (z) and po(-z) * Pd(z). 



(7) 
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FIG. 2: The dependence of the SOI amplitudes for electrons in 
the first and second subband on the DDT-parameter R, see {§). 
The initial confinement potential, at R = 0, is a square well, 
with parameters described in the caption of Fig. [3] Observe 
that «i = at R = 0.82. The curves are numerical [1^] so that 
zeros are not exactly reached. 



We proceed as follows: first we introduce a slight asym- 
metry in the positioning of the doping layers with respect 
to the center of the well, thus breaking the symmetry 
of po(z) (see Figs. [T] and [3}. On top of that, we utilize a 
technique of inverse scattering theory known as Double 
Darboux Transformation (DDT),||7|] to generate an asym- 
metric confinement potential E c (z), which satisfies the 
first of inequalities 0. 

In a DDT one changes the single particle potential 
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E c and the wave functions qff — > cp n in such a 
way that the eigenvalues e„ are preserved. This, in turn, 
offers a way to significantly modify the strengths of the 
SOI amplitudes without grossly altering the spectrum 
of the quantum well. Naturally, the invariance of the 
eigenvalues under the DDT is rigorous only if electron- 
electron interactions are neglected. However, we have 
found that, even in the presence of interactions, the vari- 
ation of the eigenvalues remains relatively small, about 
10%. 

The precise form of the transformation is[0] 



E c = Ef - 2 [<p 2 A]' , <p a = - AI 2 i, <p 2 = 



- ™(°) 



rA 
r 2 -!' 



A = 



(r 2 -l)<P? 
l + (r 2 -l)I 22 



•J -co 



(8) 



where E^ ' is the original (symmetric) confinement, and 
cp^ are the corresponding wave functions (WFs). The 
parameter R = log r measures the strength of the sym- 
metry breaking of Ef'; R = corresponds to the identity 
transformation, i.e., E c = E^ and <p„ = cpf . We as- 
sumed in © that (pi i2 (—oo) - +0 and J^ oo ip 2 (z)dz = 1. 

We further assume that Ef\z) represents a rectangular 
well parameterized by ITo (height) and a (width). 

In Fig. |2]we show the dependence of cc\ and a 2 on the 
parameter R. We see that a\ = when R = 0.82. At this 
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FIG. 3: Confinement in the conduction and valence bands at 
R = 0.82, see Fig. [2] The vertical walls correspond to the 
borders of the initial square well of width a = 400A and depth 
U = lOmeV. The left and right spacers have widths 5 L = 
0.1a and 5r = 0.04a, respectively. Both doping layers have 
widths w = 0.06a. The two lowest subbands are at energies 
£i = 4.5meV and e 2 = 8.8meV from the bottom of the original 
well. El 



value of R we have a 2 - 0.18 in the chosen units. The 
form of the confinement potential at this value of R is 
shown in Fig.[3j together with the plot of the electron en- 
velope functions for the two subbands. Such a potential 
might be realized by MBE techniques. [14] When only the 
lowest subband in this specially engineered QW is pop- 
ulated, the electrons do not suffer spin relaxation. This 
happens in equilibrium, when the energy distribution of 
electrons everywhere in the sample is given by a Fermi 
function. When the sample is driven out of equilibrium 
(e.g., by the current or by light illumination), so that both 
bands are populated, the SOI effects are activated. 

Chan ging the population of the subbands. It is well- 
known[15] that the Rashba SOI amplitude can be con- 
trolled by a perpendicular (to the QW plane) electric field. 
Below we show that an electric field which drives a cur- 
rent in the plane of the QW, can be used as a "switch" for 
the SOI amplitude, when this amplitude is initially set 
to zero by the method described in the previous section. 

Under ordinary conditions, when a current flows in an 
electron gas, in the linear response regime, strong inelas- 
tic interactions among the electrons and between elec- 
trons and phonons establish a local equilibrium distri- 
bution, controlled by a local (position-dependent) elec- 
trochemical potential p(x): 



p(x) = (1 - x)p(0) + xp(l), 
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where x is the direction along the current, L is the sample 
length, p(0) = pi = p and p(l) = pR = p(0) + eV are the 
electrochemical potentials in the left and right contacts, 
and V is the applied voltage. Similar to p{x), also the 
subband energies and the density of states depend on 



FIG. 4: The energy distribution function, /e(x*) at a point x* 
inside a current-carrying sample (see Eq. fill ) at temperature 
T = (e 2 - ei)/30. The dashed lines represent the edges of the 
conduction subbands, and the dotted line is the local electro- 
chemical potential from Eq. (9). Because of the two-step form 
of the non-equilibrium distribution function some electrons 
have energy higher than the chemical potential fi(x*) (see the 
"dashed" area). These "hot" electrons came from the right; the 
"hottest" ones have energy ~ p R and originate from the right 
contact. Note that /e(x'} - 1 for E < p L and /e(x") - x* for 
Pl S! E S! Pr- thus, the step in the distribution function /e(x') 
occurs precisely at /e(x") = x". By applying this rule one easily 
sees that the energy distributions at the edges of the sample, 
/ E (0) and /e(1), are usual Fermi functions. 



the coordinate, 

e 1/2 (x) = e 1/2 (0) + p(x) - p(0), ei, 2 (0) = e\ i2 , 

0, £ < ei{x), 
v(x, E) = < vo, £\{x) < E < £2(x), 
2v , E > e 2 (x) , 



(10) 



where vo = m/(2nh), so that the concentration of elec- 
trons n remains independent of position. 

However, in a mesoscopic sample at small temper- 
atures and high mobility values, when the condition 
DIE 2 » t -1 is satisfied, where D is the diffusion constant 

' in ' 

and Ti n is the inelastic scattering from electron-electron 
and electron-phonon interactions, the energy relaxation 
effects are so small that an electron may pass through 
the entire sample conserving its energy. This results in a 
two-step form of an energy distribution function inside 
the sample,lll6l[l7ll see Fig.H 



/ £ (x) = (l-x)/ E (0) + x/ E (l), 



(11) 



where /e(0) and /e(1) are Fermi energy distributions in 
the left and right contacts attached to the sample, respec- 
tively. When the voltage V across the sample exceeds 
the critical value V* = (e 2 — e\)le, the energy of the "hot" 
electrons from the right contact exceeds the energy of 
the bottom of the second subband on the left end of the 
device. The fraction of such "hot" electrons is controlled 
by fTTll : it is maximal near the right contact and vanishes 
in the vicinity of the left contact. More precisely, the 
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FIG. 5: Population of the subbands in the sample, see Eq. l fl2t . 
The width of the populated zone is controlled by the applied 
voltage V. For V = 1.5(e 2 _ £1), n 2 (x) approaches its maximal 
value at x' = x'/L = 0.43. B 




FIG. 6: The coordinate dependence of the SOI amplitudes in- 
side the sample. Close to the right contact, the second subband 
is not populated (see Fig.[5}; hence a.\ = 0, and electrons do not 
experience spin-orbit scattering there. fl3ll 



x-dependent population of the second subband is (see 

Kg. m 

nijx) _ |0, when e 2 {x) > ^(1), 

v )x[[i(l) - e 2 (x)] /2, when £2(x) < , 

where n\{x) — n — n2(x), and n is the total density of elec- 
trons. In fT2l we assumed that electrons spread equally 
between the bands when their energy allows them to be 
either in first or the second subband. 

In equilibrium, the second subband would remain un- 
populated everywhere in the sample since [i(x) < e 2 {x) 
for arbitrary x. In a steady (non-equilibrium) current- 
carrying state, there is a zone inside the sample, where 
both subbands are populated. The width of this zone is 
controlled by the applied voltage, see Fig. [5] The electric 
field-induced change in the population of the subbands 
modifies the Hartree potential LTh, which in turn mod- 
ifies the values of ot\ and a 2 . In particular a\ becomes 
non-zero when the voltage V exceeds a critical value V* 
{V* = 0.8mV for ^(0) - e x = 3.5meV), with the other 
parameters being specified in the caption of Fig. [3] The 
coordinate dependence of a\ and a 2 for different voltages 
V is shown in Fig. ||6}. We conclude that the Dyakonov- 
Perel spin relaxation of electrons in the lowest (majority) 
subband can be switched on by an electric field (or a 
current) in the plane of the QW. It is essential for the 
argument that this field does not modify the envelope 
functions of the electrons. The field affects a n due to 
the subband repopulation, and the effect is large only 
because <X\ had previously been fine-tuned to be zero. 

Discussion. The method described in this Letter al- 
lows in principle to quickly suppress a spin polarization 
by applying a current. A crucial point is the feasibility of 
keeping the hot electrons hot (i.e., to avoid thermaliza- 
tion). One way to ensure this is to use short wires, low 
temperatures and high mobilities, so that the condition 
D/L 2 » t^ 1 is satisfied. Failing this, another possibility 
would be to rely on the nonlinear hot electron effect, see 
p. 120 in lll8ll , whereby, in the presence of strong inelas- 



tic scattering, the distribution at the center of the sam- 
ple is a Fermi distribution with an effective temperature 

T = [T 2 + ^(eV) 2 ] 112 . If V is sufficiently large, this will 
produce the desired population of the higher subband. 
Finally, we notice that, while our treatment has neglected 
the intrinsic Dresselhaus SOI, the latter 's action may be 
compensated by a Rashba SOI having the same ampli- 
tude, so that s x - s y becomes a conserved quantity. jl9ll So 
the proposed scheme should work even in the presence 
of Dresselhaus spin-orbit interaction. 
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